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Abstract 

A chiral-perturbative consideration of the near-threshold pp ppn^ reac- 
tion indicates that the pion-rescattering term has a substantial energy and 
momentum dependence. The existing calculations that incorporate this de- 
pendence give pion rescattering contributions significantly larger than those 
of the conventional treatment, and this enhanced rescattering term interferes 
destructively with the one-body impulse term, leading to theoretical cross 
sections that are much smaller than the observed values. However, since the 
existing calculations are based on coordinate-space representation, they in- 
volve a number of simplifying assumptions about the energy-momentum flow 
in the rescattering diagram, even though the delicate interplay between the 
one-body and two-body terms makes it desirable to avoid these kinematical as- 
sumptions. We carry out here a momentum-space calculation that retains the 
energy-momentum dependence of the vertices as predicted by chiral pertur- 
bation theory. Our improved treatment increases the rescattering amplitude 
by a factor of ^^3 over the value obtained in the r-space calculations. The 
pp — > ppn^ transition amplitude, which is now dominated by the rescattering 
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term, leads to the cross section much larger than what was reported in the 
approximate r-space calculations. Thus, the extremely small cross sections ob- 
tained in the previous chiral perturbative treatments of this reaction should 
be considered as an accidental consequence of the approximations employed 
rather than a general feature. 



PACS numbers: 13.75.Cs, 13.75.Gx, 12.39.Fe 



I. INTRODUCTION 



The high-precision measurements of the total cross sections near threshold for the 
pp ppn^ reaction have invited many theoretical investigations on this process 0- 0. 
The pion production near threshold is expected to occur via the single-nucleon process (the 
impulse or Born term), Fig. 1(a), and the s-wave pion rescattering process. Fig. 1(b). In the 
conventional treatment [|T0|, the vr-A^ vertex for the impulse term is assumed to be given by 
the Hamiltonian 

Ho = (c?-V(r-7r) - -!-{a-V,r-7r}) ^, (1) 

where qa is the axial coupling constant, and f^^ = 93 MeV is the pion decay constant. The 
first term gives p-wave pion-nucleon coupling, while the second term accounts for the nucleon 
recoil effect. The s-wave rescattering vertex in Fig. 1(b) is customarily described with the 



phenomenological Hamiltonian [10 



Til = An — —ipTz-TTi/j + AiT — "^i/jt-ttxtti/j (2) 

The coupling constants Ai and A2 determined from the experimental pion-nucleon scattering 
lengths are Ai ~ 0.005 and A2 ~ 0.05. Thus, Ai ^ A2, as expected from current algebra. The 
calculations P,p!0| based on these phenomenological vertices yield cross sections for s-wave 7r° 



production that are significantly smaller, typically by a factor of ~5, than the experimental 
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values |T|. There are, however, some dehcate aspects in the calculated cross section. First, in 
Eq.(|I]), only the second term contributes to s-wave pion production. The suppression factor 
~ niT^/m^ contained in this term drastically reduces the contribution of the impulse term. 
Fig. 1(a), and as a consequence the relative importance of the two-body rescattering process. 
Fig. 1(b), is enhanced. However, the dominant A2 term in Eq.(|^) cannot contribute to the 
ppTT^TT^ vertex in Fig. 1(b) due to its isospin structure. Thus, a phenomenological calculation 
based on Eqs.(|l|) and (H) leads to highly suppressed cross sections. 

To describe the more recent theoretical developments, it is convenient to introduce what 
we call the typical threshold (TT) kinematics. Consider Fig. 1(b) in the center of mass (CM) 
system with the initial and final interactions turned off. At threshold, (go, (f) = {fn-K, 0), 
P'w = P20 = "^N, p{ = P2 = 0, so that pio = P20 = + ko = p^ = k = -p2 

with \k\ = ^JuiT^mt^ + Needless to say, even for q = 0, the actual kinematics for 

the transition process may differ from the TT kinematics due to the initial- and final-state 
interactions. Now, for the TT kinematics we have k'^ = —171^171^^., which implies that the 
rescattering process typically probes inter- nucleon distances ~ 0.5 fm. The reaction can 
therefore be sensitive to exchange of the heavy mesons, which play an important role in 
the phenomenological meson-exchange N-N potentials. Lee and Riska ||^ pointed out that 
the shorter-range meson exchanges were indeed capable of enhancing the cross section sig- 
nificantly. Meanwhile, even though the off-shell tcN amplitudes feature in the rescattering 
process {e.g., k"^ = —m^mj^ 7^ for the TT kinematics), there is no guarantee that Tii 
of Eq.(0) describes the off-shell amplitudes adequately. Hernandez and Oset suggested 
that the s-wave amplitude enhanced for off-shell kinematics could increase the rescatter- 
ing contribution sufficiently to reproduce the experimental cross sections. However, Ref. 
used phenomenological off-shell extrapolations, the reliability of which requires further 
examination. 

Chiral perturbation theory (xPT) [0,0 serves as a consistent framework to describe 
the low-energy vrN scattering amplitudes. Taking advantage of this fact. Park et al. (to 
be referred to as PM^K) and Cohen et al. carried out the first xPT calculations for the 
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pp ppn^ reaction. The results of these two groups essentially agree with each other 
on the following major points: (1) the pion rescattering term in a xPT treatment is signifi- 
cantly larger than in the conventional treatment; (2) the sign of the rescattering term in a 
XPT treatment is opposite to that obtained in the conventional approach; (3) the enhanced 
rescattering term in a xPT treatment almost cancels the impulse term, leading to theoret- 
ical cross sections much smaller than the observed values. Thus a systematic treatment of 
the off-shell vrA^ scattering amplitudes indeed drastically changes the pp — > ppir^ cross sec- 
tion. We note, however, that the calculations in Refs. which rely on coordinate space 
representation, involve potentially problematic approximations on the kinematical variables 
appearing in the energy-dependent vrA^ scattering amplitudes. Namely, in deriving the r- 
space representation of the two-body transition operator [see Fig. 1(b)], one limits oneself 
to the Feynman amplitude corresponding to the TT kinematics and Fourier-transforms it 
with respect to pi, p2, pi and P2-, while keeping all the other kinematical variables fixed at 
their TT kinematics values. Although this type of simplification of kinematics is commonly 
used in nuclear physics to derive effective r-space operators, it is expected to be much less 
reliable for the threshold pp — * ppvr^ reaction. The reason is two- fold: First, the energy- 
momentum exchange due to the initial and final-state interactions is vitally important for 
this process. Secondly, the destructive interference between the one-body and two-body 
terms implies that even a small change in the two-body term can influence the cross sections 
significantly.^ In this article we carry out a xPT calculation for the pp — ^ ppir^ reaction with 
the use of momentum representation, which liberates us from the above-mentioned kinemat- 
ical simplifications. As we shall show, this improvement drastically changes the calculated 
cross section. 

This paper is organized as follows. We describe in section 11 the effective Lagrangian of 



detailed momentum-space calculation by Hanhart et al. Q, with the use of the conventional 
Hamiltonian, indicates the importance of lifting the kinematical simplifications discussed here. 
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XPT and a formalism we use to apply to processes involving two or more nucleons. 

The details of our momentum-space calculation are explained in section 111. In section IV 
we present the numerical results and comment on their salient features. The final section V 
contains discussion and concluding remarks. 



II. NUCLEAR CHIRAL PERTURBATION THEORY 

As in PM^K 0, we use the heavy-fermion formalism (HFF) of chiral perturbation theory 
|r3| . The effective Lagrangian £ch in the HFF is expanded as 

£,h = + + £(2) + .. . . (3) 

represents a term of chiral order v with v = ('^/2) — 2, where n is the number of 
fermion lines involved in a vertex, and d is the number of derivatives or powers of m^r. In 
order to produce the one-body and two-body diagrams depicted in Figs. 1(a) and 1(b), we 



minimally need terms with z/ = and 1 R]. Their explicit forms are |12 



£(0) _ iY^iyd^JJ^d^U + (t/^ + ?7 - 2)] + N{w ■ D + gj^S ■ u)N (4) 

£(1) _ _^^NiS-D,vu}N + 2cimlNNi:i{U + -2) 
2mN 

+ (C2 — ^-^)NivufN + c^Nu-uN (5) 
8m AT 

We have retained here only terms of direct relevance for our present calculation. In the 
above, U{x) is an SU(2) matrix that is non-linearly related to the pion field and that has 



standard chiral transformation properties. We use f/(x) = y 1 — [7r(x)//7r]^ -|-iT-7r(x)//^ as 
in Ref. [|12|. The large component of the heavy-fermion field is denoted by N{x); the four- 



velocity parameter t>^ is chosen to be = (1, 0, 0, 0); D^^N is the covariant derivative; S*^ is 



the covariant spin operator and = i[^^d^^ — ^df^^'^], where ^ = yU{x) The constants 
ci,C2 and C3 have been determined from phenomenology p|JT^ and their numerical values 
will be discussed later. In practical calculations, U{x) is expanded in powers of 7r{x)/ f-,, and 
only necessary lowest order terms of this expansion are kept. 



In applying xPT to nuclei [Q, one classifies multi-nucleon Feynman diagrams into ir- 
reducible and reducible diagrams. Diagrams in which every intermediate state contains at 
least one meson are categorized as irreducible diagrams, and all others are classified as re- 
ducible diagrams. To each irreducible diagram one assigns the chiral order index v defined 
by z/ = 4 — — 2C + 2L + where E^- is the number of nucleons in the Feynman 

diagram, L the number of loops, C the number of disconnected parts of the diagram, and the 
sum runs over all the vertices in the Feynman graph fl^ . It can be shown that an irreducible 



diagram of order u carries a factor {q/Ay, where g is a generic momentum characterizing 
low-energy phenomena and A ~ 1 GeV is the scale parameter of xPT. This leads to the 
general expectation that the contributions of terms with higher values of u are significantly 
suppressed in the low-energy regime. Now, the contribution of all the irreducible diagrams 
(up to a specified chiral order) is treated as an effective operator T acting on nucleonic wave 
functions. The resulting nuclear matrix elements incorporate the contributions of the re- 
ducible diagrams. This two-step procedure may be referred to as nuclear chiral perturbation 
theory. 

Applying nuclear xPT to the present case, we write the transition amplitude for the 
pp ppiT^ reaction as 

T = ($/|r|$,), (6) 

where (|$/)) is the initial (final) two-nucleon state distorted by the initial-state (final- 
state) interaction. For formal consistency, if T is calculated up to order u, the nucleon- 
nucleon interactions that generate and |$/) should also be calculated by summing 
up all irreducible two-nucleon scattering diagrams up to order u. In practice, however, 
it is common to use the phenomenological A^-A^ interactions that reproduce measured two- 
nucleon observables. This hybrid version of nuclear xPT has been applied with great success 



to electroweak transition processes |T5[, and we shall use this phenomenological version in 
this paper. 

As discussed in PM^K, the lowest-order contributions to the impulse and rescattering 



terms come from the u = —1 and u = 1 terms, respectively. With the use of Cch in E<l-(0), 
the matrix elements in momentum spaceQ of these two terms are given by 

where pj and p- {j = 1,2) denote the initial and final momenta of the j-th proton. The 
four-momentum of the exchanged pion is defined by the nucleon four-momenta at the nNN 
vertex: kj = pj —Pj, where pj = {Ep^,pj),p- = {Epi,p-) with the definition Ep = {p"^ + 
m%y^'^. The rescattering vertex function K{k,q) is calculated from Eq.(|^): 



where k = {k^, k) and q = {ug, q) represent the four-momenta of the exchanged and final 

pions, respectively. The time component of the final pion is obviously defined by Uq = 
(g^ + m2)i/2. 

In what follows, we work with T defined by 

q- _ _|_ = q-lnip _|_ ^Resc ^j^q-j 

To specify the transition operator T completely, we need the values of the low-energy coef- 
ficients Ci, C2 and C3. In Ref. [0], these parameters were determined from the experimental 
values of the pion-nucleon a term, the nucleon axial polarizability and the isospin-even 
s-wave 7tN scattering length a"*". The results are 

ci = -0.87 ±0.11 GeV~\ C2 = 3.34 ± 0.27 GeV"\ C3 = -5.25 ± 0.22 GeV-\ (11) 

As discussed in PM^K, the value of C2 + C3 can also be extracted from the known pion-nucleon 
effective range parameter fe"*" of the low energy pion-nucleon scattering amplitude: 



^Here we normalize plane- wave states according to <p' \ p>= 6{p' — p); this normalization differs 
from the one used in Ref. 0. 
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4. 1 (i + i!!^)-(!^)-(,, + ,3-^)^. (12) 

Itt V rriN/ \ Jn J orriN m^"^ 



Since C3 in Eq. (|TT]) has been deduced directly from a^^^, a quantity known with a relatively 
high precision, we may use the value of C3 in Eq. to determine C2 from the observed 
value of 6"^. This procedure yields C2 = (4.5 ± 0.7) GeV~^ 0, which is significantly larger 
than that in Eq.(|ll|). Thus, the errors quoted in Eq.(0) do not seem to reflect the entire 
range of uncertainties in the low-energy coefficients.^ It is important to examine to what 
extent the existing ambiguities in the low-energy coefficients affect the off-shell enhancement 
of the pp — s> ppTT^ reaction. In the present work, we use as the "standard" parameter set the 
central values in Eq. ( ]TT| ) and refer to it as the parameter set I. We also use the parameter 
set II, in which the value of C2 is changed into C2 = 4.5 GeV~^. Furthermore, in developing 
a certain argument later in the text, we shall be interested in the lower end of the error bars 
in the empirical Ci. For this we introduce the parameter set III, which is identical to the 
set I except the change: Ci = -0.87 GeV"^ Ci = (-0.87 - 0.11) GeV"^ = -0.98 GeV"\ 
Thus, we will consider three sets of low-energy coefficients: 

Parameter set I 

ci = -0.87 GeV-\ C2 = 3.34 GeV-\ C3 = -5.25 GeV"^ (13) 
Parameter set II 

ci = -0.87 GeV-\ C2 = 4.5 GeV-\ C3 = -5.25 GeV"^ (14) 
Parameter set III 

ci = -0.98 GeV~\ C2 = 3.34 GeV~\ C3 = -5.25 GeV"^ (15) 

Before embarking on numerical work we make two remarks. First, T'-^^^ above represents 
only the tree-diagram contribution. Fig. 1(b); loop corrections to the u = —1 impulse term 



^ In this connection, it is worth mentioning that the reliability of the empirical value of the s-wave 
ttN scattering length a"*" has recently been questioned |16|. 
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generate transition operators of order u = 1. These additional contributions generate an 
effective ttNN vertex form factor for tlie impulse term, Fig. 1(a), but according to PM^K's 
estimate, the net effect of the loop corrections after renormalization is less than 20% of the 
leading-order impulse term. We therefore neglect the loop corrections in the present work, 
even though formal consistency requires it. We will return to this question in later work. 
Secondly, the nuclear chiral counting procedure employed above is in fact best applicable 
when energy-momentum transfers to a nucleus are small, whereas the near-threshold pp — >■ 
ppTT^ reaction involves significant energy transfers go ~ ^tt- Hence we must exercise caution 
in applying Weinberg's counting rule to this reaction. In PM^K, as far as the construction 
of the transition operators was concerned, the energy-momentum transfer due to the final 
pion was ignored {i.e., the external pion was taken to be soft, ~ 0) in order to utilize 
Weinberg's original counting rule. The physical value of was used only at the stage of 
calculating the phase space integral. The approximate nature of this method is particularly 
evident for the impulse term, where the pion-emitting nucleon is off-shell by ~ m^r. This 
off-shell nucleon must interact at least once with the second nucleon before losing its off-shell 
character. It is then sensible to treat an impulse diagram accompanied by subsequent one- 
pion exchange as an irreducible diagram, even though the original Weinberg classification 
would categorize it as a reducible diagram. Cohen et al. [|| proposed a modified chiral 
counting rule that takes account of this feature. In addition, these authors argued that the 
A degree of freedom should be taken into account explicitly in xPT since the A^-A mass 
difference ~ 2771^^ is small on the chiral scale A (see also Ref. fl^). In the present work, 
which is basically of the illustrative nature, we do not address these issues but simply use 
the transition operator given in Eq. (0) with the view to concentrating on the pragmatic 
(but in our opinion very urgent) question: For a given version of xPT transition operators, 
how important is the difference between the conventional r-space calculation and a p-space 
calculation ? 
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III. NUMERICAL CALCULATIONS 



A. Momentum space calculation 

It is most convenient to carry out the calculation in the CM frame. Referring to Fig. 1, we 
denote by p and p' the relative momenta of the two protons before and after the pion emis- 
sion, respectively. In terms of p' and the momentum of the final pion q, the momentum 
of each nucleon in Fig. 1 is written as 

— * — * 

Pi = -P2 = P, Pi = P - 2' P2 = -P - 2 ■ ^ ^ 

With the normalization of plane-wave states specified earlier, the plane-wave matrix element 
of the production operator defined by Eqs. (|7|)- (p!0[) takes the following form 

<p\q\T\p> = <p\q\T^"'^\p> + <p",g|T^"^"|p>, (17) 

where 

+ [1^2] (18) 

and 

In the above, the four-momentum transfer kj (j = 1, 2) is defined by kj = {Ej^. —E^i,pj —pj) 
with E^^ = [p'j + m?j^Y/'^ . Since we are assuming that the nuclear states are described by the 
non-relativistic Schrodinger equation, the A^-A^ potential responsible for the initial- and final- 
state interactions represents only those two-nucleon diagrams that involve no energy transfer 
between the two nucleons {i.e., three-momentum transfers only). Hence the intermediate 
nucleon energies in our treatment are given by Ej;. = {pj + m^)^/^ at each pion-nucleon 
vertex. 
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The transition amplitude of the pp — > ppir'^ reaction is evaluated by taking the matrix 
element of the production operator T defined above between the initial (x*^^'*) a-^d the final 
(x''"'*) PP scattering wavefunctions. In terms of this transition matrix element, the total 
cross section is given by 

(2vr)^ 



where pi and are the asymptotic relative momenta of the initial and final pp states, 
respectively, W = 2E^- is the total energy, Vi = 2\pi\/Ef. is the asymptotic relative velocity 
of the two initial protons, and m^^ is the z-component of the spin of the jth nucleon. Eq.(pUD 
can be calculated most easily in the partial-wave representation. We therefore expand the 
transition matrix element as 

< y a I T I y^-^^ > = 

SfLfJfMf SiLiJiMi 

X E KmM) <Pf[LfSf]JfMf\%^^M\P^[US.]J^M,>. (2l) 

where the spin-angular function of the anti-symmetrized two proton state is defined by 

yil\p,m,„m,,)= E SI ^ ^'e^'^'^'^'yiMAP) 

X < ]Xns,ms,\SMs><LSMLMs\JM> , (22) 

where 5{^ls)j is the NN scattering phase shift in the eigenchannel defined by the orbital 
angular momentum L, total spin S*, and total angular momentum J. We introduce the 
reduced matrix element using the standard convention: 

<Pf[LfSf\JfMf\Ti^^^{q)\p,[USi]J,M,> = 



Mr 



I J, I J ^ 

y -Mf M, J 



<Pf[LJSf]Jf\\TlM\\P^[L^S^]Ji> • (23) 
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The substitution of Eqs. (|2T|)-(|3D into Eq.(|OD leads to 



16 Pi Jo 

E 

LiSiJiLfSfJfl-n 



(24) 



where E^^ = {{W — oOq)'^ — pf = \J — mj^ , and the maximum momentum of the 



pion is given by = y {{W — 2mNY — m1}{{yV + 2mNY ~ mlj/AW^ . Corresponding to 
the decomposition T = T^^p + 7"Resc -j^ -g^^ (p^Tj), we have 

<Pf[LfSf]Jf\\Ti^{q) \\pi[LiSi]Ji> 
=<Pj[LfSj]Jf\\rlr{q) \\P^[L.SU> + <Pf[LjSj]Jj\\T^^^^{q) \\p,[USU> , (25) 

where each of the reduced matrix elements on the right-hand side is defined similarly to 
Eq. (^. 

For the near-threshold tt production under consideration, we limit ourselves to the case 
where the produced pion is in s-wave {It, = 0) and the final pp state is in the ^Sq state. 
Accordingly, the initial pp state is restricted to be in the ^Pq state. 

Computation of the reduced matrix elements in Eq. ( ^51) involves only the radial parts of 
the scattering wavefunctions. We calculate the momentum-space radial wavefunction from 
the half-off-shell i^'-matrix 



R{LS)J,po{.P) = i ^COs((5(L5)j) 



^{P-Po) ^^K(LS)j{P,P0,W) 



(26) 



Po Pq/uin - p'^/mN_ 

where V means taking the principal-value part of the two-nucleon propagator, and po is 
the on-shell momentum defined by = 2Epfy The i^-matrix is related to the phase shift 
by pK(^LS)j{Po,Po,W) = -tan(5(i5)j) with p = 7rpo"^Af/2. By using the well-developed 
numerical methods (see, e.g., Ref. |18|), the /^-matrix (or the scattering t-matrix) can 
be calculated directly from the potential by solving the Lippman-Schwinger equation in 
momentum space. In terms of these radial functions, the reduced matrix element for the 
impulse term is written as 

12 



xp . (-g + ^f)')5{p' - p + q/2)R^Po,pXP) (27) 
while the reduced matrix element of the rescattering term is given by 

^ < ..I'^ol II 1 1 > - / / 

where k = {ko, k) = {E^— E^'^^/2,p — p' + q/2). As will be discussed in the next subsection, 
an accurate numerical calculation of the impulse term Eq. (pT]) is possible either in momentum 
space or in coordinate space. By contrast, since the rescattering term Eq.(p8D is highly 
nonlocal, its accurate numerical calculation is possible only in momentum space. 

B. Expressions in r-space 

We now relate our p-space calculation to the r-space calculations in the literature. With 
the use of the usual normalization <f\f'>= 6{f—f') for the r-space base vectors, the radial 
functions in p- and r-space are related to each other by the standard Bessel transformation 

R(LS)j,poir) = ^'^^J P^dpjLipr)R^LS)j,poip) , (29) 
and -R(LS')j,po('^) satisfies the boundary condition 

RiLS)J,po{r) ™ ^ [cos{5(LS)j)jL{por) - sm{S(LS)j)nL{por)], (30) 

where jiix) and ni{x) are the regular and irregular spherical Bessel functions. 

As for the impulse term, the momentum-space expression, Eq. (|27|), can be easily cast 
into the coordinate-space expression with the use of Eq. ( p^ ) . The result is 

'47r J(27r)32Lj^ U ^ 
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In the g — i> limit, this expression reduces to 



1 



<pA'So]\\T^Jo\\p^?Po]> 



(27r)32m^ 



X 



9Am^ 



J drr'R^s,,p,{r){-^ + ^)R^p,A^) , (32) 



and this corresponds to Eq. (54a) of PM'^K. 

As stated, the rescattering operator is highly non-local and hence its matrix element, 
Eq. (pSD , can be evaluated exactly only in momentum space. The relation between the full 
expression, Eq. (pSD, and the approximate r-space form used in the PM^K calculation is as 

follows. We first simplify Eq. ( pSD by taking the q^O limit: 

-2vri gA 



^<Pf['So]\\T^lo\\p^['Po] 



>q^0 



^(27r)32m^ U J 
r X -\- l'\ 

X |2p' + (p - xp') log R^PoAp) 



(2vr) 



where 



X = 



2pp' 



and K is defined [cf. Eq.(|^)] as 



K 



J TT 



2Ci - (C2 - +C3) — 

8m AT 



(33) 



(34) 



(35) 



If we furthermore freeze ko, the energy variable of the exchanged pion, at the fixed value 
ko = corresponding to the threshold pion production, then we are back with the TT 

kinematics explained in section I. For the TT kinematics, 

ko=m^/2 p^ + p'"^ + 3m^/4 



X 



2pp' 



and K becomes a constant, which was denoted by Kth in Eq. (32) of PM^K: 



_ fco=m-ir/2 171 
H ^ f^th = -J^ 

J TT 



2ci - -(C2 



9a ■ 



(36) 



(37) 



We refer (as in PM^K) to the simplifications Eqs. (0), (|36|), ( pTj) as the fixed kinematics 
approximation. In the fixed kinematics approximation, the right-hand side of Eq. ( pSf ) 
becomes sufficiently simple to be recast into an r-space expression via the transformation 
Eq . (^9l) . The resulting expression agrees with the r-space form Eq.(54b) of PM^K. 
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IV. NUMERICAL RESULTS 



To calculate the reduced matrix element <p/[^5'o] 1 1 7^^™''(g) 1 1 PipPo] >■, we can use either 
Eq.(p?D or Eg. ([3T|) ; these are completely equivalent to each other. For numerical computa- 
tions, however, we find the latter more convenient. On the other hand, the reduced matrix 
element for the rescattering term, < P/[^5'o] \ \Tj^To{q) WPifPo] is calculated with the 



use of the p-space expression, Eq.(^H]). Calculations that are based on the full expressions, 
Eq. (P7D(or Eq. (PT])) and Eq. ( PS]) are referred to as full p-space calculations. For the sake 
of comparison, we also carry out calculations with the use of the approximate expressions 
discussed in the preceding section. For given reduced matrix elements, the total cross section 
is obtained via Eqs. (H), (p5|). 

The radial functions for the pp scattering states that appear in Eqs. (|28|) , (|3T|) are to 
be generated with the use of realistic nucleon-nucleon interactions. We adopt here as our 
standard choice the Argonne V18 potential [|l^. For the purpose of comparison, we also 



carry out supplementary calculations with the Reid soft-core potential pO[. As regards the 
low-energy coefficients, Ci, C2 and C3, our standard choice is the parameter set I [Eq. (p!3|)], 
and we use it throughout unless otherwise stated. 

We note here that the accuracy of the numerical integration for the rescattering term, 
Eq. (pSD, can be achieved only when the high momentum component of pp wavefunction is 
accounted for correctly. This is due to the fact that the C2 and C3 terms in n{k,q), Eq.(|), 
grow linearly with the momentum for higher momenta. We have found that a very large 
number of mesh points are needed to achieve sufficient numerical accuracies in evaluating 
the rescattering term. 

Figs. 2-4 show the total cross sections calculated with the Argonne VI 8 potential and 
with the parameter set I. It is informative to first consider the behavior of the individual 
contributions of the impulse and rescattering terms and then proceed to discuss the behavior 
of their coherent sum. Figs. 2 and 3 give the individual contributions of the impulse and 
rescattering terms, respectively, while Fig. 4 presents the cross sections due to the coherent 
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sum of these two terms. 

The cross sections in Fig. 2 have been obtained by retaining only the impulse term in the 
decomposition Eq. (^). The solid curve corresponds to the full calculation based on Eq. 
(pT|), while the broken line has been obtained with the g — > approximation, Eq. (p^). Our 
numerical results for the latter case coincide with those in PM^K. As expected, the q = 
approximation employed in PM'^K becomes less accurate as the collision energy increases. 
Although the difference between the two results is not extremely large as far as the isolated 
contribution of the impulse term is concerned, yet we should remember that, depending on 
the pattern of interference between the impulse and rescattering terms, even a rather minor 
change in the impulse contribution may affect the cross section significantly. It is therefore 
recommended to avoid the q = approximation in computing the impulse term. 

We next discuss the contribution of the rescattering term. Fig. 3 shows the cross sections 
obtained by retaining only the rescattering term in Eq. p5|). The solid curve corresponds 
to the full p-space calculation. The dash-dotted curve gives the result obtained with the 
use of the approximation Eq. (P^), in which q is fixed to be zero while the energy-transfer 
ko in the exchanged pion propagator and in the coefficient K,{k, q) is treated exactly. The 
dashed curve represents the results of the fixed kinematics approximation, viz., g = and 
fco = vtiT^jl] our numerical results for this case reproduce those obtained in PM^K. The 
difference between the dash-dotted and dashed lines clearly indicates that the non-local 
effects due to the /co-dependence enhance the cross sections by a very large factor. The 
smaller but still substantial difference between the solid and dash-dotted curve implies that 
the g-dependence of the pion propagator and of the rescattering coupling coefficient k cannot 
be neglected in a quantitative calculation of the pp — > pp-n^ cross section. The solid curve and 
the dashed curve differ by a factor as large as ~ 10. Thus the magnitude of the rescattering 
amplitude obtained in the full p-space calculation is ~ 3 times larger than that obtained in 
the fixed kinematics approximation. In Refs. the impulse and rescattering terms gave 

transition amplitudes of about the same magnitude, and this feature is also visible in Figs. 2 
and 3 (compare the dashed lines therein). In the full p-space calculation, however, it is 
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the rescattering term that is dominant, and this feature drastically changes the interference 
pattern of the two terms, as discussed immediately below. 

In all the cases shown in Figs. 2 and 3, our numerical results indicate that the reduced 
matrix element for the impulse term <pf[^So] \ \ Tj^^^{q) \ \pi[^Po] > and that for the rescat- 
tering term < Pf[^So] \ \ Ti^^^^q) ||pjpPo] > have opposite signs so that they interfere with 
each other destructively. The cross sections given by the coherent sum of these two terms 
are given in Fig. 4. The solid curve represents cxfuii, the cross section obtained in the full 
p-space calculation. The dash-dotted curve corresponds to the approximation Eq. (^). The 
cross section (Xfix.km. corresponding to the fixed kinematics approximation, if plotted in the 
same scale, would be hardly visible. The impulse and rescattering contributions in this case 
cancel each other almost perfectly, leading to an extremely small value of cr gx. km.- This is 
consistent with the previous r-space calculations f^JH- To illustrate the calculated o"fix. kin., 
we multiply it by a factor of 30 to obain the dashed curve in Fig. 4. When the approximation 
ko = is removed, the scattering term is substantially enhanced (see Fig. 2), and as a 

consequence there is only a partial cancellation between the impulse and rescattering terms. 
The corresponding cross section (dash-dotted curve) therefore becomes significantly larger 
than (Tfix. kin.- This tendency is even more prominent as we go to the full p-space calcula- 
tion by removing the g — approximation; the solid line representing cTfuu exhibits a huge 
enhancement over cxfix. kin.- 

The results in Fig. 4 indicate that the extremely small cross sections reported in Refs. 
I^J^ are an accidental consequence of the kinematical approximations used there rather than 
a general feature of calculations based on nuclear xPT. It is to be noted, however, that 
cTfuu in the present calculation is still significantly smaller than the observed cross sections. 
Our calculation here does not include the Coulomb repulsion effect in the initial and final pp 
states. The inclusion of the Coulomb effect is expected to further reduce the cross section 
0], worsening the discrepancy between the calculated and observed cross sections. 

We have already mentioned that the results presented in Figs. 2-4 are obtained with the 
Argonne V18 potential and with the parameter set I, Eq. ([T3D . To provide some measure of 
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the sensitivity of the cross section to input A^A^ potentials, we have repeated a full p-space 
calculation with the use of the Reid soft-core potential |2y], while retaining the parameter 
set I for the low-energy coefficients. The resulting cross section is given in Fig. 5 (dashed 
curve) along with the result for the Argonne V18 potential (solid curve). The difference 
between these two results is not significant. 

To examine the sensitivity to different choices of the low-energy coefficients, we have 
carried out full p-space calculations, with the Argonne VIS potential, for the parameter set 
II, Eq. (0), and for the parameter set III, Eq. ([T5|) . The dashed curve in Fig. 6 represents the 
cross section corresponding to the parameter set II. This curve is almost indistinguishable 
from the solid curve obtained with the parameter set I even though these two sets have 
significantly different values of c^- At first sight, this feature is puzzling because the large 
difference between cxfuu and cTfix. kin. seen in Fig. 4 seems to suggest that the cross sections 
should also be sensitive to the parameters contained in the vertex function K,{k,q). On 
closer examination, however, the lack of sensitivity to C2 in the full p-space calculation has 
a simple explanation. In the near-threshold region (g — ^ 0) under consideration, the C2- and 
C3- terms in the vertex function K{k, q), Eq. (^, depend on a factor = Ep — Ep> [see also 
Eq. (P5|)]. This energy factor introduces a zero in the transition operator at \p\ = \p'\, which 
drastically reduces the contributions from the C2- and cs-terms to the integration over the 
scattering wavef unctions. As a result, the rescattering term is dominated by the ci-term in 
the full p-space calculation. In contrast, in the fixed kinematics approximation{see Eq.(37)), 
the transition operator arising from the C2- and C3-terms are momentum-independent and 
of comparable magnitude but of opposite sign to the transition operator coming from the 
Ci— term. This means a substantial cancellation between these two operators at all momenta 
and hence the rescattering transition matrix element can be a sensitive function of the low 
energy coefficients including C2 in the fixed kinematics approximation calculation. This 
feature originating from the treatment of fco essentially accounts for the large difference (by 
a factor of ~ 3), discussed earlier, between the rescattering amplitudes obtained in the full 
p-space calculation and in the fixed kinematics approximation (see the discussion on Fig. 3). 
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The same feature also explains why the cross section in the full treatment is hardly sensitive 
to C2 (and C3). 

The dominance of the ci-term in the full p-space calculations implies that the cross 
section should be sensitive to an input value of ci, a coefficient linked to the pion-nucleon 
sigma term |T2[ . In Fig. 6 we compare the results of full p-space calculations corresponding 
to the parameter sets I and III, which differ only in the value of ci. The cross sections 
indeed vary substantially as we change the value of ci from the currently accepted central 
value to one end of the error bars [see Eq. (pH])]. We remark, however, that even with the 
most favorable choice of ci the calculated cross sections are significantly smaller than the 
observed values. 



V. DISCUSSION AND CONCLUSIONS 

Our calculations in this article are based on the transition operators derived from a 
particular version of the chiral effective Lagrangian, the one used in PM'^K . As mentioned 
in section II, other choices of the effective Lagrangian are possible, although the main features 
of our results may be stable against these changes. This point requires further examinations. 
A related issue is the choice of an explicit form of U (x) as a function of the pion fields 7r(x). 
Needless to say, if we calculate everything exactly, the results should be independent of 
particular choices. However, in any practical calculations one needs to expand U{x) in 
powers of 7r(x) and truncate the series, which nullifies the formal equivalence of various 
choices. Thus, in an approximate calculation such as the present one, the pp ppn^ 
transition amplitude does depend on the off-shell vrA^ amplitude, which varies when different 
forms for U{x) are used. We have employed here a particular form for U{x), and the stability 
of our results against other choices is yet to be studied. 

In our calculation, the rescattering transition matrix element corresponding to the C2- 
and C3-terms are found to be rather sensitive to the higher momentum components of the 
nuclear (pp) wave function; as a matter of fact, the relevant range of momentum is not much 
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smaller than the chiral scale A. This uncomfortable feature is in fact shared also by other 
known applications of nuclear xPT, and a satisfactory solution to this problem seems to 
require the study of terms with higher chiral orders than considered here. 

Keeping in mind all the caveats stated above, we still consider it almost certain that, in 
any reasonably realistic xPT calculations, the rescattering term dominates over the impulse 
term and their signs are opposite to each other. This implies that the heavy- meson exchange 
contributions considered in Ref. cannot be invoked as a possible mechanism to enhance 
cTfuu to bring it closer to the observed cross sections. Since it is established that these heavy- 
meson contributions have the same sign as the impulse term, the addition of these extra 
contributions to the transition amplitude obtained in the full p-space calculation would 
result in a destructive interference. Thus the heavy-meson contributions such as considered 
in Ref. suppresses the cross section instead of enhancing it. Most recently, van Kolck, 
Miller and Riska ||^ considered a yet another diagram involving p — uo exchanges. Since 
the sign of this additional contribution is unknown, one way to bring our CTfun closer to the 
experimental cross sections is to include this new contribution assigning to it the favorable 
sign. Needless to say, from a xPT point of view, this is a highly ad hoc prescription, but it 
is a possibility. 

As mentioned in section II, we have not considered here the loop corrections to the 
impulse term, corrections that lead to the form factor effect for the impulse vertex. Qual- 
itatively speaking, the inclusion of this effect is expected to reduce the contribution of the 
impulse term, further enhancing the dominance of the rescattering term. This aspect awaits 
further study. 

We now summarize. For the pp — ppir^ reaction near threshold, we have calculated the 
cross sections, using the effective transition operators derived from chiral perturbation theory 
(xPT) and employing the momentum- space representation. Our p-space calculation is free 
from the various kinematical approximations that went into the previous xPT calculations 
based on the r-space representation. In what we call full p-space calculations, we 
have retained all the energy-momentum dependence in the transition vertices as dictated 
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by xPT. The results of the previous r-space calculations are recovered by introducing into 
the present full p-space calculation the fixed kinematics approximation, which consists in 
freezing the momentum q of the final pion and the energy variable ko of the exchanged pion 
at their values corresponding to the threshold kinematics. The results of the full p-space 
calculation indicate that the fixed kinematics approximation is not justified. Specifically, 
we have found that in the full p-space calculation the contribution of the rescattering term 
[Fig. 1(b)] is enhanced by a factor of ~ 3 over the value obtained in the approximate r- 
space calculations As a result, the rescattering term now dominates over the impulse 

(Born) term, removing the near complete cancellation between these two terms that existed 
in the previous calculations. This means that the extremely small cross section reported in 
Ref. becomes much larger in the full p-space calculation [see Fig. 4]. The enhancement, 
however, is not sufficient to explain the observed cross sections. As our work is based on a 
particular version of the chiral Lagrangian used in Ref. [0], it remains to be seen to what 
extent the use of other chiral Lagrangians would affect the results. It is also interesting to 
study whether the remaining discrepancy between the observed and calculated cross sections 
can be explained in terms of diagrams involving heavy meson exchanges [Q . Thus, the near- 
threshold pp ppn^ reaction invites further detailed investigations. The principal lesson 
we draw from the present calculation is that in any serious xPT calculations one must use 
what is called here the full p-space calculation, avoiding the conventional approximations 
employed in r-space calculations.. 
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FIGURES 

FIG. 1. Diagrams contributing to pp ppir^: impulse term (a) and rescattering term (b). 
Besides the four-momenta specified in the figures, we use in the text the following notations: pi = 
relative momentum between pi and pi; pf = relative momentum between p'^ and P2! P = relative 
momentum between pi and P2- P = relative momentum between p'j^ and Tp^] W = piQ + P20 = 
total energy. 

FIG. 2. Total cross sections calculated with the impulse term alone. Solid curve: full p-space 
calculation; dashed curve: q = approximation. 

FIG. 3. Total cross sections calculated with the rescattering amplitude alone. Solid curve: 
full p-space calculation; dash-dotted curve: q = approximation; dashes curve: fixed kinematics 
approximation. 

FIG. 4. Total cross sections that include both impulse and rescattering amplitudes. Solid curve: 
full p-space calculation; dash-dotted curve: q = approximation; dashed curve: a^x. kin. ^ 30. The 
experimental points are from Ref. ^ (open diamonds) and Ref. |^ (solid diamonds). 

FIG. 5. Total cross sections obtained in full p-space calculations with the Argonne V18 potential 
(solid curve), and with the Reid soft-core potential (dashed curve). The experimental points are 
from Ref. [Q] (open diamonds) and Ref. |^ (solid diamonds). 

FIG. 6. Total cross sections obtained in full p-space calculations with the parameter set I (solid 
curve), the parameter set II (dashed curve) and the parameter set III (dash-dotted curve). The 
experimental points are from Ref. ^ (open diamonds) and Ref. |Q] (solid diamonds). 
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